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ABSTRACT 



Context. In 1985, at the end of the active phase 1977-1986, a broad (4000 kms^^) Lya line appeared in the symbiotic system CH 

Cygni that had never been observed previously. 

Aims. In this work we investigate the origin of this anomalous broad Lyo' line. 

Methods. We suggest a new interpretation of the broad Lya based on the theory of charge transfer reactions between ambient hydrogen 

atoms and post-shock protons at a strong shock front. 

Results. We have found that the broad Lya line originated from the blast wave created by the outburst, while the contemporary optical 

and UV lines arose from the nebula downstream of the expanding shock in the colliding wind scenario. 
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1. Introduction 

Symbiotic stars (SSs) were introduced as spectroscopically pe- 
culiar objects by Merril (1919). Nowadays they are interpreted 
as binary systems composed of a compact hot star, generally 
a white dwarf (WD), and a cool giant, whose mutual interac- 
tion via accretion processes is at the origin of complex spectra 
recorded at all wavelengths, from radio to X-rays. 

The first models attempted to explain the symbiotic spectra 
in the light of a binary star configuration through the nebular 
emission arising from the photoionized portion of the cool giant 
wind. The photoionizing source was played by the WD alone. 
This simple scenario was fairly able to explain the first optical 
and radio data, but became inadequate for the detailed modeling 
of line emission (Nussbaumer & Vogel 1990) and for observa- 
tions in new spectral ranges (Li & Leahy 1997). The observa- 
tional evidence that actually both stars lose mass through im- 
portant winds (Nussbaumer et al 1995) led to the colliding wind 
(CW) models (Girard & Willson 1987, Li et al. 1998, Kenny & 
Taylor 2005). A complex network of shock waves throughout the 
symbiotic systems could finally explain the observations, partic- 
ularly in the high energy spectral ranges. The coupled effect of 
shocks and photoionization was therefore accounted for in the 
modeling of the line and continuum spectra (e.g. Contini 1997). 
Moreover, the network of gas and dust nebulae within specific 
SSs could be distinguished by composite models (Angeloni at 
al. 2007 a,b,c). 

CH Cygni (CH Cyg) is amongst the most studied yet still dis- 
cussed SSs. In a recent paper we have attempted a comprehen- 
sive and self-consistent modeling of continuum and line spec- 
tra in different epochs within the CW framework (Contini et al 
2007, hereafter Paper I). 
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In particular, in 1977, CH Cyg underwent a powerful out- 
burst that lasted until 1986. Towards the end, bipolar radio and 
optical jets appeared (Solf 1987) and spectral variations were 
observed in the UV and optical ranges. A remarkable charac- 
teristic of the 1985 spectrum was the appearance of a broad, 
strong Lya emission line (Fig. 1 ), never evident in previous spec- 
tra (Selvelli & Hack 1985). This very peculiar spectral feature 
was, however, disregarded in Paper I, because it could not be ex- 
plained by the CW model alone. The overall model of CH Cyg 
has been amplified and is presented in this paper. 

We discuss the origin of the Lyff line appearance in 1985 
and suggest an alternative interpretation accounting for shock 
waves that accompany the WD outburst during the active phase. 
In Sect. 2 we briefly review the observational trend and the previ- 
ous theories suggested to explain the composite Lya line profile. 
In Sect. 3 we describe the colliding-wind scenario and the result- 
ing picture at the end of the 1977-1986 active phase (Sect. 3.1). 
Consequently, adopting Chevalier's (1982) theory on the inter- 
action of the outburst with the circumstellar medium, we apply 
the results by Heng & Sunyaev (2007, hereafter HS07) about the 
Lya broad line formation at high-velocity shock fronts (Sect. 
3.2). A discussion appears in Sect. 4, while concluding remarks 
follow in Sect. 5. 



2. The broad Lya line: observations and former 
interpretations 

During the outburst started in 1977, no Lya emission from CH 
Cyg was recorded. Then, towards the end of the active phase 
(January 1985), a strong and broad Lya emission line, as wide 
as 20 A, appeared (Figs. 1 and 2), with a strong redward-shifted 
wing, suggesting expansion velocities up to 4000 kms^^. This 
outstanding behaviour has not been repeated during the sub- 
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Fig. 1. lUE spectra of CH Cyg, showing the appearance of a 
strong and unusually broad Lye in 1985, at the end of the 1977- 
86 active phase (adapted from: MSH88, top and central panels - 
Selvelli & Hack 1985, bottom panel) 



sequent evolution throughout the active and quiescent phases. 
For instance, during the 1990-91 and 1995 quiescent stages, the 
Lya virtually disappeared, together with the hot UV and op- 
tical continuum, while during the active phase 1992-1995 the 
Lya line was similar in profile (Fig. 2), but 2-3 times weaker 
than in 1985 (Skopal et al 1998). 

The standard Lya profile has been decomposed into two - 
stable and variable - components (Skopal et al. 1998). While the 
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Fig. 2. Evolution of the Lyo- profile in the high-resolution lUE 
spectra, from 1980 to 1995. On the left, the observation date are 
indicated; the g marks the position of the narrow Lya geocoronal 
emission (adapted from Skopal et al. 1998). 



former emission has usually been attributed either to the atmo- 
sphere of the red giant (Selvelli & Hack 1985) or to some cir- 
cumbinary material (Skopal et al. 1998), the latter, i.e. the vari- 
able part of the profile that extends up to a high velocity (~ 2200 
km s'^), has been explained as a fast outflow in the vicinity of the 
hot star An asymmetric high- velocity outflow (2000 km i"' ) was 
also indicated by the broad Balmer lines occurring in the spec- 
trum on short time scales (lijima et al. 1994). Such velocities 
and variabilities are similar to those characteristic of the shock 
between the stars created by collision of the winds (Paper I). 
Therefore, they are less anomalous than the broad (4000 km s"') 
Lya line observed at the end of the 1977-1986 active phase. 

A mechanism for explaining this broadening was proposed 
by MSH88 who claimed that, according to Johansson & Jordan 
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(1984), the Lya line with an enhanced red wing, typical of line 
formation in an acceleration outflow, may be broadened by scat- 
tering by a high opacity. However, a stellar origin is difficult 
to reconcile with the sudden appearance of the broad Lya line. 
MSH88 then proposed a formation region displaced from the 
orbital plane, connected with outflowing material, while the 
B aimer lines might originate in an accretion disk. 

3. Theoretical scenarios 

We would like to explain the exceptionally broad Lya line emis- 
sion in the frame of the shock-front network in CH Cyg. The re- 
sults of the wind-collision model (Paper 1) led to a detailed phys- 
ical and morphological picture of the emitting nebulae within the 
SS. In particular it was demonstrated that the spectra depend on 
the system phase. In January 1985, at the end of the active phase 
that started in 1977, the UV and optical line ratios revealed that 
we were facing the circumstellar medium of the WD opposite the 
red giant. In this region, the dynamical consequences of the WD 
outburst can be compared with those of a supernova explosion, 
though on a different scale. Following Chevalier (1982), 8 years 
after the burst, the expanding blast wave reached a relative large 
radius and the velocity of the shock front is in the range of those 
characteristic of broad Lya line emission by charge-transfer re- 
actions (HS07). 

3.1. Results of the colliding-wind model 

Collision of the winds (Girard & Willson 1987, Kenny & Taylor 
2005) from the two-component stars leads to two main shocks: 
the head-on shock between the stars, facing the WD, and the 
head-on-back shock expanding from the system outwards (men- 
tioned in Paper I as the reverse and expanding shocks, respec- 
tively). 

According to the colliding wind model, during the active 
phase, the broad Hfi lines with FWOM of 400-1200 km i ' were 
emitted downstream of the head-on shock between the stars 
(Angeloni et al. 2007a). Although this shock is actually a stand- 
ing shock, it may be accelerated throughout the decreasing den- 
sity of the atmosphere by the massive wind from the red gi- 
ant. The broad Balmer lines emitted downstream are particu- 
larly strong because the intensity of permitted lines depends on 
the temperature of the star T, and on the ionization coefficient 
U. T, can reach more than 100,000 K during the outburst and 
U is relatively high (> 1) because the emitting gas is close to 
the hot star The broad lines decline both in intensity and width 
during quiescence because T, becomes < 30,000 K and the 
high velocity wind slows down. In a shock dominated regime, 
Rayleigh-Taylor (R-T), Richtmyer-Meshkov (R-M), and Kelvin- 
Helmholtz (K-H) instabilities at the shock fronts lead to frag- 
mentation. Adopting a filling factor between 0.001 and 1 (Paper 
I), a maximum geometrical thickness of the filaments ~ 10''' cm, 
and an average velocity of 600 km i"' during outburst, the vari- 
ability time scale is between < 1 hour and < 1 month. The time 
scale increases during quiescence depending on the velocity de- 
crease. 

The sudden appearance of the broad (4000 km ,? ') Lya line 
looks quite anomalous. Such a broad line could not even come 
from the nebula downstream of the expanding shock, which cor- 
responds to Vs< 150 A:mi '(Paper I). Interestingly, the broad 
Lya line appeared at about the same time as the line spectrum 
observed after November 1984, when the outburst was almost 
over. The spectrum revealed [OIII] 4363 and [OIII] 5007 lines, 
as well as high ionization UV emission lines, e.g. NV, CIV, SilV, 



Hell, OIII] (SelvelU & Hack 1985, MSH88). Forbidden opti- 
cal lines cannot survive the densities of 10**-10^ cm^^ that are 
characteristic of the shock between the stars. They must have 
been emitted from the nebula downstream of the shock expand- 
ing outward throughout the external region of the system, which 
in accordance with the colliding wind model (Girard & Willson 
1987), is located on the side of the WD circumstellar region op- 
posite the red giant. Consequently, the Lya line, appearing at 
this very epoch, was emitted from high-velocity gas outflowing 
throughout the extended hemispherical region on the side of the 
WD less disturbed by the dynamical effects of the red giant wind. 

3.2. The blast wave from the outburst 

We adopt MSH88 suggestion that the broad Lya emission region 
is connected with the outflow material. 

Following Chevalier (1982), we consider the interaction of 
the outburst with circumstellar matter on the assumption that is 
built up by a steady wind. If the ambient density is described by 
p oc r~* (p= 1 .4 mn n, where n is the density in number cm^^ 
and m// the mass of the H atom), the steady wind corresponds 
to s=2. The interaction of the freely expanding matter with the 
surrounding medium gives rise to a high-energy density region 
bounded by shock waves. Two shock fronts develop, one pro- 
ceeding inward in the high-density region, the other expanding 
outward in the circumstellar medium. 

The case of uniform expansion gas is described by s=2 and 
7=5/3. For s=2 the radius of the outer shock is given by the 
Primakoff' solutions: Rbw^O E/2 n A)^'^ f'^ (Chevalier 1982, 
Eq. 5), where E is the total energy, po = A R"^ (po= 1-4 m// no, 
where no is the pre-shock density of the gas upstream), and t is 
the time elapsed from the burst. 

This equation is valid for times longer than the time of 
change t, = 0.677 M^^^/A E'^^, between that of the interaction 
of freely expanding matter with the surrounding medium and the 
following one, i.e. when the flow tends toward the self-similar 
solution for a point explosion in a power-low density profile 
(Sedov 1959). Here M stands for the ejected mass. 

The velocity Vs- dR/dt is constrained by the observed 
FWOM of the Lya line 8 years after the burst. After some al- 
gebra we obtain no=l-45 lO"''^ E. The mass ejected by the WD 
in the CH Cyg system during the 1977-86 outburst is about a few 
(2-3) 10"'' Mq (Taylor & Seaquist 1985). The total energy is half 
thermal and half kinetic (Chevalier 1982). The high velocity ob- 
served after January 1985 presumes that the velocity of the ejecta 
was ~ 14,000 km i^'at t=l year If all the ejecta had a velocity 
of 14,000 km i"', the associated kinetic energy would amount to 
6 lO''^ erg, the total energy E=1.2 lO'"' erg, and consequently no 
- 0.2 cm^^. This in turn gives RBvt'=2.25 10'^ cm. 

To check this result we calculate for instance, no close to the 
WD (hoito) ^t a radius Ro,^^ - 2 lO'-' cm. We choose the distance 
that was found for the standing shock between the stars (Paper 
I) by modeling the continuum SED observed in May 1985. 
Although the interbinary region shows mixing of the winds from 
the stars, the standing shock facing the WD reflects the composi- 
tion of the WD (e.g. Contini 1997). From mq™ =«0s„ Rw/Rq,™ 
where noen, - 0.2 cm^^, we find no,™ = 2.5 10^ cm^^ in good 
agreement with the preshock density no -5 10^ cm'^ adopted 
to explain the continuum SED at that epoch (Paper I, Table 1). 
Constraining t, within a maximum radius of ~ 2. 10'^ cm and 
adopting no-Q.2 cm^^, we find tj= 5.2 years, confirming that 
Chevalier's model is valid in this case. 
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We can now adopt the theory of HS07, who treated the blast 
wave shocks in the SNR case. Their results show that broad 
Balmer and Lyman lines are produced by charge transfer reac- 
tions between the post-shock protons and the ambient atoms. A 
population of post-shock atoms follows with a broad velocity 
distribution (broad neutrals). For Vs> 500 km i^'the broad neu- 
trals can produce Lya that is blue- or red- shifted by resonance 
with the stationary atoms, hence providing an escaping way for 
the protons. For shocks with Vs> 4000 kms^\ the luminosity 
ratio rj^ya/Ha is > 10 (HS07 Fig. 1), so the Lya line would be 
strong. 



Interestingly, the energies calculated from the temperature (~ 
2.94 lO'' K) downstream of this strong shock front correspond to 
~ 250 keV, not far from the gamma-ray energy range, while en- 
ergies of ~ 21 keV (2.4 10** K) correspond to y,=4000 km s'K 
We suggest that the observations of symbiotic star outbursts be- 
yond hard X-rays might lead to interesting results. 

The escape velocity of the WD is v^sc =(2 G Mwd/Rivd)'^^ 
~ 1.6 lO'^/R,;,^^, adopting a mass Mwd ~ IM© ■ We recover a 
velocity of 14,000 A:mi"' similar to what is predicted at early 
times for the Lya line FWOM profile, adopting a WD radius Rwd 
~ 1.36 10** cm in the CH Cygni system. 



4. Discussion 

As previously mentioned (Sect. 3.1), the broad Lya appeared at 
about the same phase as the expanding shock within the collid- 
ing wind scenario. The radius of the expanding nebula Kexp ~8 
10^^ cm results from modeling the observed spectral energy dis- 
tribution of the continuum (Paper I). Notice that the radius of 
the blast wave Kbw calculated in Sect. 3.2, is about three times 
that of Kexp, while the physical parameters, yj=4000 km .s^'and 
no=0.2 cm^^ of the blast wave, and ^,=150 km .s^'and «o=10^ 
cm^^ of the expanding shock, are very different. This is not sur- 
prising, considering that the blast wave stems from the outburst, 
while the expanding shock derives from the wind collision. The 
two shock fronts will hardly interfere even in the orbital plane, 
because R-T and K-H instabilities lead to fragmentation of mat- 
ter at the shock fronts with filling factors < 0.01 (Paper I; Contini 
& Formiggini 2001). 

A velocity of ~ 4000 km .s^'was evident from the Lya line 
profile after January 1985 in CH Cyg. The other UV and optical 
lines showed narrow FWHM profiles. In contrast, UV lines in- 
dicating expansions of 3000-4000 km s"^ were observed after the 
1985 eruption of the recurrent nova RS Ophiuci (Snijders 1987; 
Shore et al. 1996) and also optical lines after the last outburst 
in 2006 (Bode et al. 2007 and references therein) showed high- 
velocity expansion. The UV and optical lines strong enough to 
be observed can be emitted downstream of a strong shock if the 
preshock density is high enough to speed up the cooling rate (oc 
n^) downstream. Consequently the temperature will drop below 
10* K, leading to recombination. However, the densities are con- 
strained by the critical densities for collisional deexcitation of 
the different ions, which are relatively low (< 10*" cm^^) for for- 
bidden lines. In order to emit strong enough forbidden lines (e.g. 
[OIII]), a shock with ^5=4000 km ,? 'should propagate through- 
out a medium with no ^ 10"* cm^^ The lines are emitted from 
the gas beyond the temperature drop characteristic of shock- 
dominated regimes downstream, at a distance from the shock 
front of > 7 10'* cm. Lower densities correspond to stronger 
forbidden lines at larger distances from the shock front. Strong 
[FeVII]6087, [FeX] 6375, and even [FeXI] 6986 are predicted 
at such high V,. The UV lines are generally permitted or semi- 
forbidden, therefore they can also be emitted at higher densities. 
The OVI 1034 and CIV 1500 lines will be particularly strong. 
These predictions imply a preshock magnetic field of Bo =10"^ 
Gauss. 

Such densities are not unusual in the circumstellar medium, 
as they were found by modeling the [Oil] and [Nil] lines dur- 
ing the active phase 1998-2001 of CH Cyg (Paper I). The high 
preshock density would indicate that the shock is interacting 
with mass lost by either the progenitor star or by the SS in a 
previous burst. However, such high velocities were never seen in 
the line profiles of CH Cyg, except for the broad Lya in 1985. 



5. Concluding remarlts 

In 1985, at the end of the active phase 1977-1986, a broad (4000 
kms'^) Lya line was observed that had never been present in 
previous spectra from the symbiotic system CH Cygni. 

We have noticed that the broad Lya line appeared contem- 
poraneously with the optical-UV spectrum emitted downstream 
of the shock created by collision of the winds from the stars. 
This shock expands throughout the extended circumbinary re- 
gion located on the side of the WD circumstellar region opposite 
to the red giant; as a result, the Lya line originates somewhere in 
the hemispherical region opposite the WD, where the dynamical 
consequences of the burst are less affected by the red giant wind. 
The overall situation is similar to that of an SN explosion. This 
suggests that the broad Lya emission line and the other optical 
and UV lines observed at the same phase could be emitted from 
different shocked nebulae. We consider that the broad Lya line 
stems from the WD outburst. 

Applying the theory developed by Chevalier (1982) for Type 
II supernovae to the interaction of the WD outburst with circum- 
stellar matter, we have found that the expanding blast wave had 
reached a radius Rbw ~ 2.25 lO'^ cm 8 years after the burst for 
a shock velocity of 4000 km s^\ We then applied the theory de- 
veloped by HS07 for high-velocity shock fronts in SNR, namely, 
the broad Lya line is produced by charge transfer reactions be- 
tween the blast wave post-shock protons and the ambient pre- 
shock atoms. For shocks with Vs> 4000 km s^\ the luminosity 
ratio Liya/Ha is > 10 (HS07), so the observed Lya line is strong. 

The energy involved with the outburst is E ~ 1.2 10^* erg, 
and the ambient density consistent with a velocity of 4000 
kms'^is no ~ 0.2 cm'^. Higher velocities of about 14,000 
fern ,? 'predicted by the Sedov solution at early times, lead to 
temperatures in the downstream region of such a strong-shock 
front, corresponding to emission in the near gamma-ray fre- 
quency range. Such emission was not observed at that time, first 
because of technical inadequacy and also because, during the 
active phase, the WD circumstellar region opposite the red giant 
only beaome visible in 1985 when the broad Lya appeared. The 
velocities had already decreased to about 4000 km ,?"'. Actually, 
the SEDs in the VBU range during the 1977-86 period, presented 
in Paper I (Fig. 6), show that the black body flux from the hot star 
does not appear 

According to previous results (Paper I), the hard X-rays are 
emitted from a small region between the two component stars, 
while the soft X-ray are emitted from the extended circumbinary 
region. The results of this paper suggest that the hardest radiation 
(~ 250 keV) comes from the WD circumstellar region close to 
the WD on the opposite side of the giant star, while hard X-rays 
could also be observed at distances < 0.08 pc from CH Cyg. 

Finally, the complex light curves of CH Cyg (Eyres et al. 
2002) show that the nebulae created by collision of the winds 
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and the dusty shells ejected from the red giant, expanding out- 
ward beyond the binary system, lead to temporary obscuration 
(Paper I). We suggest that the shock front at a relatively large 
radius corresponding to the blast wave may also contribute to 
obscuration episodes. 
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